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ABSTRACT 
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The interaction of the geomagnetic field and a neutral ionized 
corpuscular solar stream is investigated. According to the theory of 
Chapman and Ferraro, the interaction is such that the geomagnetic field 
is confined to the interior of a cavity in the stream into which the 
particles of the stream do not penetrate. An approximate solution is 
presented for the complete shape of the boundary of the geomagnetic field 
for the case where the magnetic dipole axis is normal to the free -stream 
velocity vector. The calculations were made with the aid of an IBM 7^90 
electronic computer. 
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INTRODUCTION 

It is generally believed (see, e.g.. Chapman [i960], Ferraro [1960a], 
and Dungey [1958] for recent accounts) that the earth's magnetic field 
does not extend infinitely in all directions, but is bounded on the 
sunward side as a result of its interaction with a neutral ionized stream 
of particles emanating from the sun. (Dissenting views have been summarized 
by Alfven [1961I.) Hie interaction of the neutral solar stream with the 
geomagnetic field is such, according to the theory of Chapman and Ferraro, 
that the field is confined to the interior of a cavity in the stream into 
which the particles of the stream do not penetrate. Although much of the 
earlier work is concerned with the manner in which the boundary forms and 
evolves (see Ferraro [1960a] and Chapman and Kendall [1961] for recent 
developments), the present paper is concerned with the determination of 
the ultimate shape and location of the boundary for the idealized case in 
which the earth is deeply immersed in a uniform steady stream. 

The steady-state Chapman -Ferraro problem is formulated in the following 
way. The magnetic field inside the boundary is considered to consist of 
the geomagnetic dipole field plus a component induced by electric currents 
in the boundary. It is supposed that the component of the total field 
normal to the boundary vanishes, and that the square of the total (tangential) 
field intensity at the boundary is proportional to the normal component of 
the momentum of the incident solar particles. The shape and location of 
the boundary and the total magnetic field inside it are then given by that 
solution of the magnetic field equations which satisfies the stated boundary 
conditions. The exact solution has not been given for this problem in three 



dimensions, but the corresponding problem in two dimensions has been 
solved analytically by Zhigulev and Romishevskii [ 1959 ] j Hurley [1961], 
and Dungey [1961I. 

Beard [i960] simplified the three-dimensional problem by relinquishing 
the condition that the component of the total magnetic field normal to 
the boundary vanishes, and applying instead, the approximation that the 
magnetic field at the boundary is two times the tangential component of 
the dipole field. He then used a power series technique to calculate the 
shape of the boundary of the geomagnetic field on the daytime side, for 
the case where the dipole axis is perpendicular to the free -stream velocity 
vector . 

Spreiter and Briggs [1961, 1962a, 1962b] employed a similar analysis 
and calculated traces of the boundary in the plane of the dipole axis and 
free -stream velocity vector for arbitrary inclinations of the dipole axis. 
They used the related approximation, suggested by Ferraro [1960b], that the 
total magnetic field at the boundary is 2f times the tangential component 
of the dipole field, where f is a constant. No information was supplied 
concerning the proper value for f in the three-dimensional problem. It 
was shown, however, that the distance along the sun-earth line to the 
boundary in the corresponding approximate two-dimensional problem is about 
5 percent too large if f is unity. If f is taken to be O.91, the 
coordinates of the entire approximate boundary agree well with the results 
indicated by the exact two-dimensional solutions cited in the foregoing 
paragraph . 
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The present work is an extension of that done by Beard [i960], and 
Spreiter and Briggs [1961, 1962a, 1962b]. The coordinates of the entire 
boundary of the geomagnetic field are calculated for the case in which 
the dipole axis is perpendicular to the free -stream velocity vector. 

The mathematical derivations have been given by the authors cited here 
and are not repeated in the present paper. The notation used by Spreiter 
and Briggs is retained for convenience. 

It may be observed that the dimensions of the boundary, as calculated 
here using representative stream conditions, are much smaller them have 
been inferred from the magnetometer experiments of Pioneer I and Pioneer V. 
(See Sonett, et al. [i960], and Coleman, et al. [1960I.) The dimensions of 
the boundary can be increased by substantially decreasing the velocity or 
number density of protons in the solar stream; or, as has been indicated 
by Spreiter and Briggs [1961, 1962a], by the inclusion in the theoretical 
model of a ring current as has been proposed by Smith, et al. [I960I on 
the basis of data from Explorer VI and Pioneer V. The effect of a ring 
current has not been considered in the present calculations. It is expected, 
however, in the light of results given by Spreiter and Alksne [1962], that 
the general features of the shape of the boundary would not be greatly 
affected by inclusion of a ring current , but that the dimensions of the 
boundary would be increased considerably. 

MATHEMATICAL FORMULATION OF THE PROBLEM 
Preliminary Discussion 

The determination of the shape of the boundary of the geomagnetic 
field, and the total magnetic field H Inside it, involves the solution 
of the magnetic field equations div H » 0 and curl H ■ 0 . The earth’s 
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magnetic field is represented by a three-dimensional dipole singularity 
at the origin (the center of the earth) . The normal component of H 
vanishes at the boundary, and Dungey [1958] has shown (see Spreiter and 
Briggs [1962b] for clarification of contrary views expressed in Spreiter 
and Briggs [1961a]) that the total (tangential) field at the boundary, 

H s , may be expressed mathematically by the relation 

H ^ 

= 2mnv 2 cos 2 f (l) 

8it 

The quantities m, n, and v are mass, number density, and velocity of 
the protons of the solar stream, and i|r is the angle between the free- 
stream velocity vector and an outward normal to the surface. The 
condition cos i|r < 0 must hold on the boundary. 

It is a property of the boundary value problem described above that 
the field H can vanish only at isolated points on the boundary. These 
points are designated neutral points. It follows from equation (l) that 
cos \|r vanishes, and the boundary is therefore parallel to the stream, 
at these points. 

Beard [i960] dropped the condition that the normal component of H 
vanishes at the boundary and replaced it with the approximate condition 
that H s *» 2H-fc, where is the tangential component of the geomagnetic 

dipole field Ep at the boundary. The closely related approximation, 
suggested by Ferraro [1960b], 

H s = 2fEfc (2) 

where f is a constant, was used by Spreiter and Briggs [1961, 1962a, 
1962b] and is employed in the present work. 
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The differential equation that defines, approximately , the shape of 
the boundary is obtained by substitution of equation (2) into equation (l). 
For the case where the dipole axis is normal to the free-stream velocity 
vector, it is 
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where p = r/r G and 

sVe _ Af 2 Hp Q \ 1/6 

r ° \l6rtmnv 2 / a \l6 rtmnv 2 / 

The variables r, 9, and cp are spherical coordinates that are fixed with 
respect to the geomagnetic dipole axis (see fig. l) . The quantity Mp, 
the magnetic moment of the dipole, is equal to a 3 Hp Q , where a represents 
the radius of the earth, and Hp^ is the magnitude of Hp at the magnetic 
equator on the surface of the earth. The quantities in equation (4) are 
in cgs units, and numerical values for a, Hp Q , and m are 6.37X10 8 cm, 
0.312 gauss, and 1.67*L0 -24 gram, respectively. 

The quantity r D is the geocentric distance along the sun-earth line 
to the boundary of the geomagnetic field. The effect of variations of n 
and v on r Q is shown in figure 2 for two values of the constant f. 
Representative values of v and n are 500 km/sec and 10 protons/cm 3 . 

These lead to values for r G of 7*6 earth radii for f = 1 and 6.9 earth 
radii for f = O.75. 

If attention is confined to the plane 9 = ±jt/2 (x = 0) where the 
derivative dp/d9 is zero by consideration of symmetry, equation (3) 
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reduces to an ordinary differential equation that can be solved analytically. 
The trace of the boundary in this plane is illustrated in figure 3* The 
front portion is circular, p = 1, and the upper portion is defined by 
the relation p cos 9 » (3/2 2/3 )p 3 /(l + P 3 ) • 

The problem in three dimensions is to find the solution of equation (3) 
whose trace in the plane 9 = ±it/2 is as shown in figure 3. It will be 
seen, subsequently, that the boundary in the upper half -space (z > 0) is 
composed of two intersecting surfaces. One will be called the main surface. 
It contains the circular trace p = 1, and that portion of the upper trace 
that lies in the half -plane 9 = -jt/2 (i.e., z > 0, y < 0) . The other, 
to be called the polar surface, contains the small portion of the upper 
trace that lies in the half -plane 9 = Jt/2 (i.e., z > 0, y > 0) . This 
segment of the upper trace is indicated by broken lines in figure 3. The 
coordinates of the boundary of the geomagnetic field are then determined 
by solving the initial value problem defined by equation (3) together with 
the traces of the solution in the plane 9 = n/2. 

Equation (3) may be put into a form that renders it more amenable to 
standard numerical procedures by solving for the derivative bp/ciq>; that is. 



The functions F x , F 2 , and F 3 are 
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F, * sin 0 - S .05-g. 

P be 


F 2 = -L fsin 0 + 2_£os_e dfi) 

p 3 \ p be/ 
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The numerical integration of equation (5) is discussed in the next section. 


NUMERICAL INTEGRATION 


The coordinates of the two portions of the boundary of the geomagnetic 
field are determined by forward integration from the traces in the plane 
9 “ Jt/2 toward larger values of cp. The variable p is extrapolated over 
intervals Acp in surfaces of constant 9 by use of an iterative extrap- 
olation process based on Euler’s method for solving ordinary differential 
equations. (See, e .g . , Kunz [1957]-) Equation (5) is used to calculate 
the derivatives dp/ 69 and at each new value of cp, bp/be is evaluated by 
numerical differentiation. 

The initial step away from the plane 9 = jt/2 is effected by means of 
the truncated series 

P = P Q + 49(dp/d9) D + ( 8 2 p/ 89 2 ) o + . . . 

* p Q (l + a (7) 

which is similar to the one used by Beard [i960]. Here p Q represents 
the trace of the boundary in the plane of 9 = n/2. The constant a - a.(o) , 
which is proportional to the second derivative of p with respect to 9, 
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is evaluated by substituting equation (7) into equation (5). It develops 
that a is a solution of the quadratic equation 

A f iSo 1 f * 

0 J + 2^ cos0 J a + -r = o 

(8) 

where S 0 = (8 p/80 )q >=1t / 2 > and the functions f f z , and F 3 are simply 

and F 3 (eqs. (6)) evaluated at cp = n/2. In the case of the main 
surface, where p Q 85 1 and 5 0 = 0, the desired solution of equation (8) is 

a = i j^(2 + 3 sin 2 ©) - J(2 + 3 sin 2 ©) 2 - k sin 2 ©^ (9) 

After a is computed by means of equation (8), the values for p at 
9 = n/2 + A9 are readily calculated by means of equation ( 7) • The deriva 
tive 8 p/89 here is simply 2ap Q and 8p/8© may now be calculated by 
numerical differentiation. Thus the trace of the boundary, as well as its 
derivatives, is completely determined in the plane 9 = n/2 + Acp. 

The integration process for succeeding intervals proceeds in the 
following fashion: The values for p and 8p/89 at station 9^ are 

inserted in the linear extrapolation formula 

P i+1 = Pi +Ap(8p/89) ± (10) 

in order to get a first estimate for p at station 9i +1 = 9i + ^9. The 
derivative (8p/8©)i +1 is computed by numerical differentiation, and then 
equation (5) is used to evaluate the derivative (8p/89)^ +i . Finally a 
refined value for Pi +1 is determined by means of the iterative 
extrapolation formula 
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p i+i = P 1 + (^P/2)^(Sp/a9) i + (dp/dcp) i+x (ll) 

The derivatives (Sp/^0)^ +1 and (dp/dcp)^ +1 are now recomputed as before , and 
the process is repeated from equation (10) to integrate over the next 
increment Acp. 

The forward integration technique outlined here has been used to find 
the coordinates of traces of the boundary of the geomagnetic field in 
planes cp = constant in the range Jt/2 < 9 < 3 it/ 2 and 0 < 0 < n/2. This 
range of values of 9 and 0 is sufficient by reasons of symmetry to 
establish the entire boundary. In any plane 9 = constant the traces of 
the two surfaces intersect, and the boundary is defined by the exterior 
portion of the intersecting traces. 

RESULTS AMD DISCUSSION 

The coordinates of the main and polar surfaces have been calculated 
at 5° intervals of 9 in the range 90° < 9 < 270°. The increment and 
range of variation of 0 were 5° and 0 < 0 < l40° in the case of the main 
surface , and 2° and 0 < 0 < 3O 0 in the case of the polar surface . It is 
clear that the main surface possesses mirror symmetry in the equatorial 
plane, and it will be seen subsequently that the range of interest of 0 
in the polar surface is 0 < 0 < 20°. The calculations were carried out 
beyond the range of interest of the variable 0 for both surfaces, how- 
ever, in order to minimize the effects of end-point inaccuracy inherent 
in the procedure used to evaluate the derivative 3p/c)0 in planes of 
constant 9. It may be observed that the results do not change significantly 
with the use of increments somewhat larger or smaller than those employed 


here . 
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It develops that the numerical technique described in the foregoing 
section is not adequate for calculating the coordinates of the main sur- 
face near the plane 9 = 270°. The derivative dp/69 becomes excessively 
large in this region, with the chosen mesh size, and therefore the extrapo- 
lation process is not sufficiently accurate. The procedure works well up 
to 9 = 255 ° t however, and these results along with the known analytic 
solution in the plane 9 - 270° (see fig. 3)> are stiff icient to determine 
the entire boundary. 

The numerical results are presented graphically in figures and 5 • 
Figure ^ shows traces of the boundary of the geomagnetic field in the 
planes 9 = 90° > 105°, 120°, . . ., 270°. The trace for the plane 9 = 270° 
is the analytic solution as given in figure 3* Figure 5 is an isometric 
view of the boundary. The curve designated S in this figure represents 
the intersection of the main and polar surfaces, and the point labeled N 
corresponds to the neutral point in the exact solution. 

It is of interest to examine the deterioration of the results with the 
introduction of certain simplifications in the integration procedure. It 
is found, for example, that the power series starting process may be omitted 
with no major effect on the results. A comparison is shown in figure 6(a) 
of calculations made with and without the starting series. The differences 
are seen to be confined essentially to traces of the main surface. If, on 
the other hand, equation (ll) is deleted, then the extrapolation process is 
purely linear. A comparison is shown in figure 6(b) of calculations made 
with the full iterative process, and the linear process with and without the 
starting series. It is seen that if the starting series is retained, linear 
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extrapolation produces results that are in substantial agreement with the 
results obtained using the iterative process. If, on the other hand, the 
starting series is dropped, the linear process yields irregular results 
that would probably be regarded as unsatisfactory for most purposes . 

The present results may be compared with those given by Beard [ i960 ] . 
Traces of the boundary in planes 90 ° < 9 < 180 ° are reproduced in figure 7 > 
and Beard's results are indicated by broken lines. It is seen that the 
two sets of computations agree well except in a region near the polar axis. 
These differences result from corresponding differences in the trace of 
the boundary in the meridian plane cp = 9O 0 that have been discussed pre- 
viously by Spreiter and Briggs [1961, 1962a, 1962 bl. Beard has informed 
the present authors privately that the latter differences result from the 
fact that the present authors base their analysis on the consistent use of 
the approximation given by equation (2), whereas he uses another approxima- 
tion in the vicinity of the polar axis. Beard, and also Spreiter and Briggs 
[1961], have in addition made refined calculations for the trace in the 
equatorial plane by solving the ordinary differential equation that results 
when 9 is set equal to 90 ° and 6p/d0 to zero in equation ( 5 ). These 
results are in almost perfect agreement with the present calculations. 

As noted in the introduction, Spreiter and Briggs [1961, 1962a, 1962b] 
showed that the coordinates of the boundary given by the exact and approximate 
solutions of the two-dimensional Chapman-Ferraro problem can be brought into 
close agreement if the value for the constant f is taken to be 0 . 91 . The 
principal remaining difference is that the discontinuous slope in the vicinity 
of the neutral point in the approximate solution appears as a smoothly turning 
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curve in the exact solution. The exact solution has not been given for 
the three-dimensional problem, but according to an estimate of Beard [ i960] 
the geocentric distance on the sunward side to the approximate boundary is 
too great by about 7 "to 11 percent . A major portion of this discrepancy 
can be rectified if f is chosen in the range 0.7 to 0 . 8 . There remains, 
however, the discrepancy associated with the discontinuity of slope along 
the curve of intersection of the main and polar surfaces. 

In the exact solution of the Chapman -Ferraro problem, the field 
strength is small in the vicinity of the neutral point, and vanishes at 
the neutral point. It follows from equation (l) that the slope of the 
boundary with respect to the free-stream direction is small in this region, 
also. These remarks apply to the approximate solution as well, except that 
the field strength does not vanish at the point that corresponds to the 
neutral point. Although the percentage errors may be large, the absolute 
magnitude of the errors in slope should be small or moderate. These errors 
would tend, in any case, to be of secondary importance since the coordinates 
of the boundary are calculated by a process involving the integration of 
surface slopes. It is anticipated, therefore, that the discontinuous slope 
along the curve S represents a local failure of the approximate method 
of solution of the problem. It is concluded, however, that with the appro- 
priate choice for the value of the constant f, the coordinates of the 
approximate boundary should be in substantial agreement with the coordinates 
that would be given by the exact solution of the Chapman -Ferraro problem. 
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Fig. 1. 
Fig. 2. 

Fig- 3- 
Fig. 4. 

Fig. 5- 
Fig. 6 . 

Fig. 7 . 


The coordinate systems. 

Values of r Q for various velocities and number densities of the 

solar stream. 

Analytic solutions in the plane 9 = ±jt/2. 

Traces of the boundary of the geomagnetic field in planes 
of constant cp. 

Isometric view of the boundary of the geomagnetic field. 

Comparisons of results obtained by iterative and linear extrapolation 
processes with and without the starting series. 

Comparison with Beard's results. 
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